HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

等热流密度加热垂直下降液膜的研究

利用等热流密度加热条件下降膜流动的三维模型方程进行线性稳定性分析和数值模拟。线性稳定性分析表明，模型方程在小到中等Reynolds数下都适用，并且流向不稳定性增长率随着Reynolds数和Marangoni数增加而增加，展向不稳定性增长率则随着Marangoni数增加而增加，随着Reynolds数增加而减小，流向和展向对扰动波数都存在一个不稳定区间。三维数值模拟表明，在等热流密度加热条件下，液膜在随机扰动的情况下最终会形成带孤立波的三维溪流状结构，液膜与气体的换热也因溪流状结构的出现而加强；在随机扰动的基础上引入占优势地位的展向最不稳定扰动会使得换热增强，液膜会提前破裂；在随机扰动的基础上引入占优势地位的流向最不稳定扰动时，液膜的换热会增强，但不会提前破裂；在随机扰动的基础上同时引入占优势地位的流向和展向最不稳定扰动时，换热会加强且液膜会提前破裂。     

一种考虑界面相物性影响的界面单元法

界面是由复杂的界面相简化而成的,界面破坏实际是界面相材料的破坏。数值计算为了方便,如经典模型和内聚力模型等,都把很薄的界面相作无厚度化处理。导致只能考虑界面的面力,而无法考虑界面相内的应力(平行于界面方向的应力)。使界面失效准则先天性地排除了界面相内部应力的影响,从界面相材料失效机理的角度来看这是不够严谨的。本文将界面相材料等效为一种弹性连续体,由界面本构关系推导得到了一种新的界面单元。该单元具有界面参数易确定、对界面相物性可以进行等效描述等优点。通过商用有限元软件ABAQUS和用户子程序UEL实现了数值分析,并与直接物理模型的数值模拟结果进行对比,证明了本方法的简便及准确性。通过对不同界面相厚度结构的进一步分析,探讨了本文方法的可行范围。     

Molecular insight into structures of monocationic and dicationic ionic liquids in mica slits

Structures of ionic liquids (ILs) 1-decyl-3-methylimidazolium bis(trifluoromethanesulfonyl)azanide ([C₁₀mim][TFSA]) and 1-decyl-dimethylimidazolium bis(trifluoromethanesulfonyl)azanide ([C₁₀(mim)₂](TFSA)₂) in different-sized mica slits have been investigated using molecular dynamics simulations. Ion density and angular distributions for monocationic IL [C₁₀mim][TFSA] were analysed to elucidate the IL structures under different surface charges and especially their changes in the direction perpendicular to the surfaces. [C₁₀mim][TFSA] formes in bilayers, compatible with existing models of ILs with long alkyl chains. For dicationic IL [C₁₀(mim)₂](TFSA)₂, cations adjacent to the mica surface tend to stay parallel to the surface with both positively charged rings absorbed. While near the centre of the slit, dications show the weak tendency of orientation distribution, more random than [C₁₀mim]⁺ ions. Structures of [C₁₀(mim)₂](TFSA)₂ cannot be described by bilayer models. Additionally, the in-plane arrangement of [C₁₀mim][TFSA] is more ordered when K⁺ ions completely neutralise the negative charge of the mica surface, and [C₁₀mim]⁺ ions tend to be located in hexagonal mica lattices with two aluminium atoms in replacement of silicon atoms. [TFSA]^? ions are constrained by the neighbouring K⁺ ions absorbed onto mica lattices.     

Back to the future: A retrospective assessment of model‐based scenarios for the management of the shrimp fishery in French Guiana facing global change

While the number of models dedicated to predicting the consequences of alternative resource management strategies has increased, instances in which authors look back at past predictions to learn from discrepancies between these and observed developments are scarce. In the past decades, the French Guiana shrimp fishery has experienced shrimp market globalization and decreasing levels of shrimp recruitment due to environmental changes. In 2006, a bio‐economic model of this fishery was developed to simulate its possible responses to economic and environmental scenarios up to 2016. Here, we compare here these predictions to the observed trajectories. While the number of active vessels corresponds to that which was predicted, the estimated shrimp stock does not. Important driving factors had not been anticipated, including a general strike, natural disasters, and the end of the global financial crisis. These results show the importance of participative approaches involving stakeholders in the co‐construction and shared representation of scenarios. Recommendations for resource managers

• Effective fisheries resources management and a fortiori, the capacity of the fisheries to adapt to global change, requires understanding of both ecological and economics dynamics.
• The temporal trajectory of the trawling shrimp fisheries has been well monitored, and the decline of both stock and fleet is understood regarding ecological and economic changes: Changes in the environmental conditions of shrimp recruitment, and oil price increase and selling price decrease.
• However, our bio‐economic modeling work showed that, even with a good understanding of the dynamics explaining past trajectories, unpredictable events (strike, natural disasters…) have acted as other key driving factors altering the capacity of the model to represent possible futures.
• These results led us to recommend a better integration of the expertise of social and political scientists in developing models of bio‐economic systems to increase the quality of scenario predictions, and to argue for more participative approaches involving the stakeholders.

     

基于有限体积虚拟区域方法的两相流动直接模拟

为提高计算效率，提出有限体积法离散下的虚拟区域颗粒两相流动直接模拟方法.在控制方程中加入相应的虚拟区域源项，保证了颗粒内部的刚体运动特性.该源项中含有颗粒信息部分及流体信息部分.在每次迭代后，对源项中的流体信息部分进行更新，从而更好地保证颗粒内速度的刚体分布.计算静止颗粒圆柱绕流及单个颗粒的沉降过程，验证了算法的准确性.     

一类正交空间填充设计的构造

空间填充设计在计算机试验中应用十分广泛,当拟合回归模型时,正交的空间填充设计保证了因子效应估计的独立性.基于广义正交设计,文章给出了构造二阶正交拉丁超立方体设计和列正交设计的方法,新构造的设计不仅满足任意两列之间相互正交,还能保证每一列与任一列元素平方组成的列以及任两列元素相乘组成的列都正交.当某些正交的空间填充设计不存在时,具有较小相关系数的近似正交设计可作为替代设计使用.设计构造的灵活性为计算机试验在实践中的广泛应用提供了必要的支持.     

基于LBM方法的风力发电机尾流结构仿真

针对传统CFD数值计算方法难以实现风力机动态旋转及其旋转状态下的流固耦合计算,本文结合格子玻尔兹曼(LBM)方法易于处理动态复杂边界的特点及大涡模拟(LES)方法在非稳态涡流结构捕捉上的优势,采用LBM-LES联合方法进行三维风力发电机整机气动性能及尾流结构仿真研究,同时采用尺度自适应方法对尾涡结构进行跟踪和精细化计算。针对NREL PhaseⅥ型试验机进行模拟,得到了与实验结果吻合的流动形态及尾流结构演变规律,分析了尾流区速度演变规律并对比了不同亚格子湍流模型对计算结果的影响.